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We report on an extensive characterization of the cracking noise produced by charcoal samples 
when dampened with ethanol. We argue that the evaporation of ethanol causes transient and ir¬ 
regularly distributed internal stresses that promote the fragmentation of the samples and mimic 
some situations found in mining processes. The results show that, in general, the most fundamen¬ 
tal seismic laws ruling earthquakes (Gutenberg-Richter law, unified scaling law for the recurrence 
times, Omori’s law, productivity law and Bath’s law) hold under the conditions of the experiment. 

Some discrepancies were also identified (a smaller exponent in Gutenberg-Richter law, a stationary 
behavior in the aftershock rates for long times and a double power-law relationship in productivity 
law) and related to the different loading condition. Our results thus corroborate to elucidate the 
parallel between seismic laws and fracture experiments caused by a more complex loading condition 
that also occurs in natural and induced seismicity (such as long-term fluid injection and gas-rock 


outbursts in mining processes). 

PACS numbers: 89.20.-a, 89.75.Da, 62.20.mt, 05.40.-a 

Earthquakes and fractures of materials are phenomena 
deeply connected under the crackling noise idea mm , in 
which systems under slow perturbation respond through 
discrete events with a huge variety of sizes. The most 
fundamental seismic laws also emerge in laboratory-scale 
experiments related to the fracture of materials [M] 
and have been recently reproduced by numerical discrete 
element simulations of porous materials [IZIIIH]. In these 
experiments, an external and constant loading is applied 
to the material and the system’ response is usually ob¬ 
tained by recording acoustic emissions. A constant and 
compressive loading is considered the most suitable anal¬ 
ogy with natural seismicity, since the main stresses un¬ 
derlying tectonic earthquakes are considered compressive 
and stationary m- In fact, a very complete parallel be¬ 
tween the acoustic emissions produced by a porous ma¬ 
terial under constant (uniaxial) compression and earth¬ 
quakes was recently reported by Baro et al. [501. How¬ 
ever, there exist other important situations related to 
natural and induced seismicity that do not fit the pre¬ 
vious conditions. This is the case of seismic events pro¬ 
duced by long-term fluid injection |2TJ [22] and gas-rock 
outbursts caused by the release of gas that is common 
and represent a serious threat in coal mining [531 HI] ■ 

In these situations (where the loading is internal, tran¬ 
sient and irregular), a complete parallel between the 
cracking noise of materials and the fundamental seismic 
laws has not been established yet, despite the consider¬ 
able interest in mining processes. Here we design a simple 
experiment that captures the previous features. Specif¬ 
ically, we study the acoustic emissions of charcoal sam¬ 
ples dampened with ethanol. At room temperature, we 
observe that the ethanol is absorbed through the pores 
of the samples and soon evaporates, creating different 
and irregularly distributed internal stresses that promote 
the fragmentation of samples and somehow mimic the 


situations found in mining. We show that these acous¬ 
tic events fulfill the Gutenberg-Richter law [55H55] (with 
a power-law exponent smaller than those reported for 
earthquakes) and the unified scaling law for the recur¬ 
rence times between events |29I4fl4j (with parameters very 
close to those reported for small mine-induced seismic¬ 
ity [55]). We also characterize the sequence of after¬ 
shocks/foreshocks, where the Omori’s decay [55Ij37| is 
observed to hold only for short times (~6 seconds), from 
which these rates display a stationary behavior. Still on 
the aftershock sequences, we investigate the productivity 
law |38j , where a double power-law relationship between 
the number of aftershocks and the energy of the trigger¬ 
ing mainshock is found (the first power-law exponent is 
much smaller than those reported for earthquakes, while 
the second is in the range of earthquakes). We also find 
that the relative difference in energy magnitude between 
the mainshock and its largest aftershock approaches the 
value of 1.2 as the mainshock energy increases (that is, 
an approximate quantitative agreement with the Bath’s 
law [39j ). Thus, in general, we verify that the fundamen¬ 
tal seismic laws hold in a fracture experiment caused by 
an internal, non-stationary and irregular loading; how¬ 
ever, our results also reveal some significant differences (a 
smaller exponent in Gutenberg-Richter law, a stationary 
behavior in the aftershock rates for long times and a dou¬ 
ble power-law relationship in productivity law), which we 
attribute to the different loading condition. 

In the experiment, charcoal samples (Fig. for do¬ 
mestic use (~200 g), made of Eucalyptus sp., are damp¬ 
ened with -30 ml of ethanol (for domestic use, hydrated 
with 7% water), and most of it is absorbed through 
pores of the samples. The results we report are based 
on six samples of about the same size. After -- 5 min¬ 
utes, the samples start to produce the cracking noise that 
is recorded by a condenser microphone (Shure Microflex 
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MX202W/N) positioned at ~20 cm from the samples with 
a sampling rate of 48 kHz. The samples emit sound for 
~30 minutes and during the processes they also crack, 
ending-up very fragmented (Fig. Figure shows 

the normalized sound amplitudes (that is, the original 
amplitudes divided by the aggregated standard devia¬ 
tion), A{t), recorded during this process. Notice that 
the sound emissions occur in discrete events with differ¬ 
ent magnitudes. The normalized energy associated to an 
event i is evaluated (analogously to other fracture experi¬ 
ments) via Ei = where I{t) = 

tinii (tendi) represents the start (end) time of the event. 
The value of tini; is chosen as the time for which I(t) ini¬ 
tially exceeds a threshold /min, whereas tend; represents 
the time for which I{t) stays below Jmin for more than At 
seconds. All results presented here were obtained with 
/min = and At = 0.1; however, different values for 
these parameters (we have tested for /min from 10~^ to 
6 X 10"^ and At from 0.025 to 0.5) do not change our 
results. The location time associated with an event i is 
dehned as U = (tend; - linii)/2. We have verified that the 
rate of activity r(t) (number of events per minute) dis¬ 
plays an approximate power-law decay in the beginning 
of the process followed by a nearly stationary behavior 
(Fig. [^). In our analysis, we have dropped out the 5% 
initial and final events in order to keep the activity rates 
nearly stationary (Fig. SI). 



Time, t (minutes) 




FIG. 1. Schematic description of the experiment. Fig¬ 
ures (a) and (b) show pictures of a sample (Sample ^2) 
immediately before and at the end of the acoustic emis¬ 
sions (see a video of this process and a selection of large 
events at http://youtu.be/8Y_uJWVeU7w and http://youtu. 
be/ILCMgFtQjnc ). In (b), we note several fissures caused by 
the cracking process. Figure (c) shows the normalized sound 
amplitudes A{t) recorded during this process. Figure (d) 
shows that the rate of activity r{t) (number of events per 
minute) displays an initial power-law decay (r(t) ~ f-o.6±o.i ^ 
for t < 1 min) followed by a nearly stationary behavior with 
average (r) = 100 ± 5 (see also Fig. SI). 


We start by evaluating the probability distribution for 
the energies E. Figure shows this distribution for one 
of the samples, where it exhibits a remarkable power- 


law behavior compatible with the Gutenberg-Richter law, 
that is, P{E) -- l/A^, over several decades. Figure]^ 
shows the values of /3 estimated via maximum likelihood 
method for different low energy cutoffs E*. We note that 
/3 is quite stable over E*. In order to assign a charac¬ 
teristic exponent 0) to each sample, we have evaluated 
the average of (3 over E*. The six samples yield values 
for P in the range [1.27-1.33], which are smaller than the 
(3 « 1.67 observed for earthquakes PSIEH] and in the Barb 
et al. m experiment (/3 rs 1.40). 




FIG. 2. The Gutenberg-Richter law. Figure (a) shows the 
probability distribution of the energies E using data from one 
experiment (Sample #2). The dashed line represents a power- 
law decay where P{E) ~ 1/E^ with f) = 1.30. Figure (b) 
shows the values of the power-law exponents j3 obtained via 
maximum likelihood method as a function of a lower energy 
cutoff E* , that is, considering only events with E > E* . The 
error bars are 95% bootstrap confidence intervals. The value 
of P is the weighted average of P over E* , where the weights 
are chosen to be inversely proportional to the lengths of the 
confidence intervals. The results for other samples are very 
similar (Figs. S2 and S3). 

Another important aspect of earthquakes is related to 
the time intervals r between events above a lower bound 
energy Amin (also called recurrence or waiting times). 
Bak et al. [2^1 have proposed that after accounting for 
the spatial location of the events, the distributions of r 
collapse into a single curve. Corral [301132] has argued 
that the occurrence of earthquakes differs from region to 
region and proposed an extension to the Bak et al. pro¬ 
cedure, by including the local rates of seismic activities 
r^y in the scaling operation. Thus, the distributions of 
T become self-similar, P{t) = rxyf{rxyT), where f{x) is 
a scaling function. When r^y is time-dependent, 
exhibits different power-law regimes that are almost uni¬ 
versal across several different seismic regions; whereas 
for Txy nearly stationary, f{x) is usually adjusted by 
a gamma distribution, f{x) oc x'^~^ exp{-xlb), where 
7 and b are fitting parameters. In fact, as proposed 
by Saichev and Sornette |33j . this behavior is an emer¬ 
gent property of aftershock superposition that holds in 
real seismicity under certain conditions [31]. Figure 
shows the distributions of r obtained from one of the 
samples and by considering several values of Amin, where 
it is clear that A(r) depends on Amin- Figure]^ shows 
the same distributions (for all samples) rescaled by the 
mean rates of activity (r). We observe a good col- 
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lapse of the distributions and that the gamma distri¬ 
bution is a reasonable fit to the average behavior with 
7 = 0.69 ± 0.08 and b = 1.50 ± 0.12. These values are very 
close to those reported for earthquakes (7 = 0.67 ± 0.05 
and b = 1.58 ± 0.15 [JT]) and small mine-induced seismic¬ 
ity (7 = 0-74 ± 0.02 and b = 1.35 ± 0.06 [35]) — see also 
Refs. mdOHIS- 



FIG. 3. Self-similarity of the recurrence times and the univer¬ 
sal scaling law. Figure (a) shows the probability distributions 
of recurrence times r with E > i?min using data from one ex¬ 
periment (Sample #1). Each curve is associated with a value 
of Emin, as indicated by the color code. Figure (b) shows 
the distributions rescaled by the mean rates of activity (r) 
(gray lines) using data from all experiments (Samples #1-6). 
The black circles are the window average over all distribu¬ 
tions and the error bars are 95% bootstrap confidence inter¬ 
vals. The solid red line is the gamma distribution adjusted 
to the average distribution via ordinary least squares method 
(the parameters are shown in the plot). Similar results are 
observed for each sample (Figs. S4 and S5). 

We now focus on quantifying the Omori’s law in our 
data. The Omori’s law [35H37] establishes that the num¬ 
ber of aftershocks per unit of time (^ms) decays as a 
power-law function of the elapsed time since the main- 
shock tins, that is, i?a(ims) ~ l/^ms' T^e value of p for 
earthquakes differs from one catalog to another (proba¬ 
bly due to different tectonic conditions), usually lying in 
the range [0.9-1.5] [35]; its value also depends on the 
magnitude of the mainshock m- In fracture experi¬ 
ments, Hirata [3] showed (for basalt) that the value of 
p decreases during the fracturing process and Bard et 
al. [3D] reported a quite stable Omori’s decay (about six 
decades) with p = 0.75 ± 0.10. In our case, we define 
the mainshock events as those with energy E^s in the 
range [IQt-lOt^^] (with j = -4, -3,..., 2) and a sequence 
of aftershocks is the events following the mainshock until 
another mainshock event is found. We calculate the af¬ 
tershock rates i?a (tms) as a function of the elapsed time 
since the mainshock, tms, averaging over all events in the 
same energy window. Figure]^ shows i?a(tms) for all en¬ 
ergy ranges and samples that we analyzed. Our results 
show that the Omori’s decay (of about two decades) with 
p = 0.87± 0.01 only holds for short times (tms ^ 0.1 min), 
from which a stationary behavior is observed for Ra{tms)- 
This stationary behavior indicates that late aftershocks 
occur randomly in time, such as in stochastic process 
with no memory. Very similar results are obtained for 


the foreshock rates (Fig. &)■ 



FIG. 4. The Omori’s law for aftershocks and foreshocks. Fig¬ 
ure (a) shows the number of aftershocks per unit of time (af¬ 
tershock rates), Ra{tms), as function of the time to main- 
shock, tms, employing data from all experiments (Samples 
#1-6). Figure (b) shows the analogue plot for the foreshocks 
(foreshock rates, Rf{tnis), versus the time before the main- 
shock, tms). The mainshocks have been defined as events with 
energy in the range lO-’-lO^^^, with j = -4,-3,...,3. Each 
gray curve is an Omori’s plot for one of the samples with 
the mainshocks in one of the energy ranges. Eigures S6 and 
S7 show the results for each sample, identifying the energy 
range of the mainshocks. In both plots, the circles (red for 
aftershocks and blue for foreshocks) are the window average 
over all curves and the error bars are 95% bootstrap confi¬ 
dence intervals. We observe power-law decays (of about two 
decades) for 7?a(tms) and i?/(tms) for tms < 0.1 min followed 
by a plateau-like behavior. The black lines are power-law 
functions, i?{oj}(tms) ~ tiS, adjusted to the average behav¬ 
iors for tms <0.1 min via ordinary least squares method. The 
power-law exponents p are shown in the plots. 

The productivity law states that the number of after¬ 
shocks Na{Ejns) triggered by a mainshock of energy 
is related to E^g via IVa(^'ms) ~ with a « 0.8 for 
earthquakes |38j . In order to quantify this law in our ex¬ 
periment, we count the number of aftershocks Na{E^s) 
that a mainshock of energy E^^s triggers. Figure]^ shows 
.^a(7?ms) versus Ejns for all aftershock sequences (defined 
as in the Omori’s analysis) in log-log scale, where (despite 
the scatter) a significant dependence is observed (Pearson 
correlation of « 0.7). This figure also shows the window 
average of these data, from which the relationship be¬ 
tween Na{Ejns) and Ej^s becomes clear: for E^^s < 10 , we 
have Na{Etns) ~ with a = 0.28 ± 0 . 01 ; whereas for 
i?nis > 10 , we find another power-law with a = 0.81 ±0.06. 
Thus, the first power-law exponent is similar to that re¬ 
ported by Bard et al. [3D] {a « 0.33), while it is smaller 
than the ones reported for earthquakes {a e [0.7-0.9] [3D] 1 
and also for creep of ice single crystals (a Ri 0.6 i)- On 
the other hand, the second power-law exponent is anal¬ 
ogous to earthquakes, but it should be considered much 
more carefully because it only accounts for about 3% of 
the data in a region where the relationship Na{E^s) ver¬ 
sus Ems displays a large scatter. 

Still on the aftershock sequences, we address the Bath’s 
law [3D], which states that the relative difference in en¬ 
ergy magnitude (that is, logE) between the mainshock 
and its largest aftershock is (on average) close to 1 . 2 , 
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regardless of the mainshock magnitude. To do so, we 
calculate the relative difference in energy magnitude be¬ 
tween a mainshock and its largest aftershock (AM = 
logi^ms - loggia, where Ei^ is the energy of the largest 
aftershock) as a function of mainshock energy E^s- Fig¬ 
ure shows the average (over all samples) of this relative 
magnitude (AM) as function of the Flms- We note that 
(AM) is systematically smaller than 1.2 for small values 
of Finis! however, AM approaches a constant plateau (for 
Eras ~ 10"^) as the mainshock energy Eras increases. This 
plateau is statistically indistinguishable from the Bath’s 
law predictions. To our knowledge, it is the first time that 
this law is studied for acoustic emissions experiments. 
The ETAS (epidemic type aftershock sequence) model 
presents a similar behavior for (AM) [53] and only for 
Eras ~ 10^ (0.8 < /3 < 1.0 in the model) that (AM) « 1.2. 
Empirical observations of the Bath’s law for earthquakes 
usually report large fluctuations for AM estimated from 
individual aftershock sequences |40j . Furthermore, aver¬ 
aged values of AM also show deviations of the Bath’s 
predictions for Eras ^10^ m, which are associated to 
lower magnitude cutoffs in earthquake catalogs — a sit¬ 
uation that cannot be ruled out in our experiment. 



Mainshock energy, E,„s 


.(b). 

- 

- 

- • 

• 

fd’l 

- 


10“® 10“* 10”® 10”® 10"’ 10° io’ 10® 

Mainshock energy, 


FIG. 5. The productivity law and the Bath’s law. The gray 
dots in (a) show the number of aftershocks NaiEms) triggered 
by a mainshock of energy Ems (in log-log scale). We have ag¬ 
gregated data from all experiments (Samples #1-6) and the 
results by samples are quite similar (Fig. S8). The red circles 
are window averages and the errors bars are 95% bootstrap 
confidence intervals. For Ems < 10, the average behavior of 
Na{Ems) is adjusted (via ordinary least squares) by a power- 
law relationship Aa(i5ms) ~ E^^, where a = 0.28 ± 0.01 (black 
solid line). For greater mainshocks (Ams > 10), the average be¬ 
havior starts to deviate from the previous power-law and can 
be adjusted by another power-law with a = 0.81 ± 0.06 (blue 
dashed line). Figure (b) shows the average value of the rela¬ 
tive difference in magnitude ((AM)) between the mainshock, 
log Ems, and its largest aftershock, log Aia, as a function of the 
mainshock energy Ems (the a:-axis is in log scale). The red 
circles are the results obtained with data from all experiments 
(individual samples show a very similar behavior, see Fig. S9) 
and the error bars are 95% bootstrap confidence intervals. 
Notice that {AM) approaches the value of 1.2 (dashed line) 
as EniB increases. 


We have thus presented an extensive characterization 
of the acoustic emissions of charcoal samples dampened 
with ethanol, aiming to establish a parallel between seis¬ 
mic laws and a fracture experiment where the loading 


(caused by absorption and evaporation of ethanol) is in¬ 
ternal, transient and irregular. We have found that the 
most fundamental seismic laws are, in general, valid in 
our experiments. However, some discrepancies with the 
case of earthquakes and fracture experiments under con¬ 
stant and external loading were also observed, nominally: 
a smaller Gutenberg-Richter exponent, a stationary be¬ 
havior in the aftershocks/foreshocks rates for long times 
and a double power-law relationship in the productivity 
law. We believe that the main cause of these discrepan¬ 
cies is the different loading condition in our experiment. 
The internal stresses are irregularly distributed across 
the samples and may create several cracking sites acting 
approximately independently. This possibility partially 
explains the observed discrepancies: simultaneous events 
yield large values of energy, which contribute for a longer 
tail in the energy distribution and for a small power- 
law exponent; cracking sites operating independently also 
corroborates with the aleatory behavior observed long 
times in the aftershocks/foreshocks rates, the Omori’s 
decay is also shorter for large values of the mainshock en¬ 
ergies (Fig. S6); finally, the crossover behavior observed 
for large values of mainshock energies in the productive 
law can result from a superposition of mainshock events. 
Another possibility is that the hydrated ethanol may pro¬ 
mote environmentally-assisted crack growth (stress cor¬ 
rosion) that also lead to acoustic emissions. We find very 
hard to direct verify these possibilities; however, stress 
corrosion usually happens in longer time scales (com¬ 
pared with our experiment), which would produce a small 
contribution for the acoustic events. We further believe 
that the simplicity of our experiment may trigger direct 
investigations related to previous discussions as well as 
the study of different solvents, sample sizes and charcoal 
materials. 
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FIG. SI. Number of events per minute (activity rate) for each sample. Notice that all activity rates are nearly stationary, 
exhibiting no clear increasing or decreasing tendency. 




FIG. S2. Probability distributions of the energies E for each sample. Each dashed line is a power-law decay with the power-law 
exponent P shown in the plots. Notice that all samples display approximately the same behavior. 
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FIG. S3. Power-law exponents /? obtained via maximum likelihood method as a function of a lower energy cutoff E* for each 
sample. The error bars are 95% bootstrap confidence intervals. The values of j3 are the weighted average (with weights inversely 
proportional to the lengths of the confidence intervals) of P for each sample. Notice that the exponents P are quite stable over 
E* for all samples. 



Recurrence times with E> E^j^, r (minutes) 


FIG. S4. Probability distributions of the recurrence times r with E > i?min for each sample. Each curve in each panel is 
associated with a value of Emin, as indicated by the color code. Notice that all samples display approximately the same 
behavior. 





































































































































Scaled recurrence times, r (r) 


FIG. S5. Probability distributions of the recurrence times rescaled by the mean rate of activity (r) for different values of ISmin 
(indicated by the color code) for each sample. In all cases, we observe a good collapse of the distributions. The black circles are 
the average over all distributions for each sample and the error bars are 95% bootstrap confidence intervals. The solid green 
lines are the gamma distribution with same parameters employed in Fig. 3b (7 = 0.69 ± 0.08 and b = 1.50 ± 0.12). 
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FIG. S 6 . Number of aftershocks per unit of time, Raitma), as function of time to the mainshock, tms, for each sample. The 
mainshocks have been dehned as events with energy in the range 10^-10^^^, with j = -4, -3,..., 3, and color code indicates the 
rage employed. In all plots, the black circles are the average over all curves of the panel and the error bars are 95% bootstrap 
confidence intervals. The black lines show the same power-law function employed in Fig. 4a, that is, Raitma) ~ t^aj with 
p = 0.87 ± 0.01. We note the power-law regime is small for large values of mainshocks energies. 
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FIG. S7. Number of foreshocks per unit of time, i?/(tms), as function of time before the mainshock, tms, for each sample. The 
mainshocks have been defined as events with energy in the range 10^-10^^^, with j = -4, ~3,..., 3, and color code indicates the 
rage employed. In all plots, the black circles are the average over all curves of the panel and the error bars are 95% bootstrap 
confidence intervals. The black lines show the same power-law function employed in Fig. 4b, that is, ~ with 

p = 0.85 ±0.01. 
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Mainshock energy, E^s 


FIG. S8. Number of aftershocks NaiEj^s) triggered by a mainshock of energy in each sample (gray dots). In each panel, 
the red circles are window average valnes and the error bars are 95% bootstrap confidence intervals. The black solid lines 
are the same power-law relationship, Na{Enis) ~ Sms, with a = 0.28 ± 0.01, shown in Fig. 5a as well as the blue dashed lines 
represent a power-law with a = 0.81 ± 0.06. 
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Mainshock energy, E„ 


FIG. S9. The red circles show the average values of the relative difference in magnitude {{AM)) between the mainshock, 
logi?ms, and its largest aftershock, logFiia, as a function of mainshock energy Ems (for each sample). The horizontal axis is 
in log scale. The error bars are 95% bootstrap confidence intervals. Notice that, for each sample, these relative differences 
approach the value of 1.2 (dashed lines) as Ems increases. 






































































































